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ABSTRACT 

We specify the range to which perturbations penetrate a planetesimal system. Such 
perturbations can originate from massive planets or from encounters with other stars. 
The latter can have an origin in the star cluster in which the planetary system was 
born, or from random encounters once the planetary system has escaped its parental 
cluster. The probability of a random encounter, either in a star cluster or in the 
Galactic field depends on the local stellar density, the velocity dispersion and the time 
spend in that environment. By adopting order of magnitude estimates we argue that 
the majority of planetary systems born in open clusters will have a Parking zone, 
in which planetesimals are affected by encounters in their parental star cluster but 
remain unperturbed after the star has left the cluster. Objects found in this range 
of semi-major axis and eccentricity preserve the memory of the encounter that last 
affected their orbits, and they can therefore be used to reconstruct this encounter. 
Planetary systems born in a denser environment, such as in a globular cluster are 
unlikely to have a Parking zone. We further argue that some planetary systems may 
have a Frozen zone, in which orbits are not affected either by the more inner massive 
planets or by external influences. Objects discovered in this zone will have preserved 
information about their formation in their orbital parameters. 

Key words: planetary systems; minor planets, asteroids: general; open clusters and 
associations: general 


1 INTRODUCTION 

Planetary systems seem to be composed of one or more stars, 
orbited by about a dozen planets and many minor bodies 
| |Galilei|[l632p . The latter can be roughly divided into hun¬ 
dreds of moons and dwarf planets, and many millions of 
planetesimals. The objects closer to the star seem to be or¬ 
ganized in a disk-like structure in which also the planets 
reside, and which flares to become spherical at larger dis¬ 
tance from the stellar host. 

This view of planetary systems is heavily based on the 
Solar System (see the review [Adams|2010p , but so far its gen¬ 
erality cannot be excluded, because each of the several thou¬ 
sands planetary systems known today ( Howard]|2013 1 seem 
to comply to this characteristic. This view is also supported 
by our limited understanding of the formation of planetary 


systems (Kokubo & Ida 2002 
et al.|2013| )! 


Bouwman et al. 2008 Ida 


The orbits of planets and minor bodies within a few stel¬ 


lar radii are affected by tidal evolution (Zahn 19771. Once 


the star leaves the main-sequence, copious mass loss starts 
to affect the entire planetary system ( [Veras et al.|[20TT| . 
We refrain from discussing these complexities here, but con¬ 
centrate on the dynamically affected regime: sufficiently far 
away and sufficiently early in its evolution to remain unaf¬ 
fected by the stellar host. 


The effect of perturbations on the Solar System either 
from the local planets, external perturbations from the birth 
cluster or even from the Galaxy have been studied quite ex¬ 
tensively (for a few recent studies see e.g. Kaib & Quinn 


|2008[ [Brasser et al.||2012[ |Schwam^|2014| and references 
therein). We start with a discussion on the Solar System 
in §|^and generalize in §|^ 


2 PERTURBING THE SOLAR SYSTEM 

2.1 The effect of internal perturbations from the 
planets 

Apart from mass loss and tidal evolution the inner regions 
of the Solar System are most strongly affected by dynam¬ 
ical interactions with the giant planets. We can calculate 
the range of dynamical reorganization caused by the widest 
massive planet in a planetary system by adopting its apoc- 
enter distance. For the Solar System, this range is currently 
determined by the planet Neptune, which affects the orbits 
of minor bodies to a distance a ^ 30 AU. Within this dis¬ 
tance the planets in the Solar System have caused major 


changes in the orbital distributions of the minor bodies (Lev- 
jison et Sj]|2008| ). We could argue that within a distance of 
30 AU/(1 — e) the minor bodies in the Solar System quickly 
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Figure 1. Fragility of the Solar System while still a member of 
its birth cluster. The shades represent the probability distribution 
rienci for fho number of encounters that perturb the Solar System 
at a given semi-major axis a and eccentricity e. For the birth 
cluster we adopted the mass of 2 x 10^ Mq and virial radius of 
2 pc to remain constant over a time scale of 200 Myr. The solid 
black curve gives the semi-major axis up to which Neptune can 
perturb orbits. The blue dashed curve gives the distance to which 
star Q ( [Jilkova et al.|2015| l, that passed the Sun with an impact 
parameter of 320 AU, has perturbed the Solar System. According 
to the color scaling along the right edge of the figure, such a 
close encounter could have occurred roughly once while the Solar 
System was a member of the star cluster. The two bullet points 
give the orbital parameters of Sedna < |Brown et al.|[2004[ | and 
2012 VPii 3 | |Trujillo & Sheppard|2014[ |. 


lose memory of the mechanism that brought them in these 
orbits. Here e is the eccentricity of the planetesimal. 


2.2 The effect of external perturbations from the 
parental star cluster 

We adopt the view that all stars are born in a clustered envi¬ 
ronment ( |Lada &: LadapoM] ). The densities of these environ¬ 
ments var y enormously from ^ Istar/p c^ to more than 10 ® 
stars/pc^ ( Portegies Zwart et al.|[2010 |. The time that the 
Solar System with minor bodies remains in the environment 
determines the degree by which it is dynamically affected 
by stellar encounters (see the review by Davies et al. ] 200 ^ . 
Those encounters play a major role in the evolution of these 
environments. The degree by which the planetary system is 
affected by such encounters depends on the duration and the 
intensity of the exposure. Further internal reorganization of 
the perturbed planetary system enables the external pertur¬ 
bations to propagate, generally on a much longer time scale, 
to the inner parts of the planetary system. 

The way in which a planetary system is affected by dy¬ 
namical encounters depends on the mass of the encountering 
star, its impact parameter with respect to the host star, the 
velocity v and the direction with respect of the planetesimal 
disk ( Steinhausen fc Pfalzner|2014 K This complicated com¬ 
bination of parameters and their mutual relations in terms 
of the degree of perturbations for the planetary system can 
be summarized in a cross section (cr). 

Most important for preserving the integrity of the plan¬ 
etary system is its eccentricity. Here we adopt a lower limit 


to the perturbation of the eccentricity of Se — O.le to be 
fatal for the particular planetesimal. The cross section of ec¬ 
centricity perturbation of a planetary system around a star 
with mass m, through an encounter with another star of 
mass M has been estimated by means of integrating small- 
N systems and averaging over the various encounter angles 
I Li fc Adams|20l5 l: 

(rr) = (1 /b) (rTsingle) T f b {^binary) ■ ( 1 ) 


Here (Tsingie is the cross section for encountering a single star 
and Ubinary is for binaries. The parameter fb is the binary 
fraction, which is between 0 and 1. We adopt fb — 0.5. In the 
adiabatic regime where the encounter velocity is compara¬ 
ble to the orbital velocity and which is suitable for encounter 
in star clusters, both cross sections have a similar form ( [lI| 
fc Adarns||2015| (with subscript X=single or X=binary de¬ 
pending on the configuration of the encountering object): 

(ifo^) ®-P(Kl-ef)).(2) 


Here v is the relative velocity for which we adopted the 
velocity dispersion in the cluster and a is the orbital semi¬ 
major axis of the planetary system before the encounter. 
The post encounter eccentricity, er = e -I- 5e. Eq.[^was cal¬ 
ibrated for initially circular orbits, but we apply them here 
also for eccentric orbits. According to Li fc Adams] ( |2015 1 
the cross sections do not depend much on the pre-encounter 
eccentricity (but see |Heggie fc Rasio|1996p . 

The parameters no, b and 7 depend on the binarity of 
the encountering object. For a single star ctq — 1000 AU^, 
6 = 8/5 and 7 = 6/5, whereas for a binary ag ~ 4050 AU^, 
6 = 4/3 and 7 = 7/5 (see Li fc Adams|2015 |. With this cross 
section the local stellar density n and the velocity dispersion 
(v) we can calculate the encounter rate: 


F = n{(j){v). 


(3) 


In Figurej^we present the expected value for the num¬ 
ber of encounters in the Sun’s parental star cluster. Here we 
adopted the cluster parameters derived by [Portegies Zwart | 


(20091: a total mass of about M^i = 2 x 10 '” Mq, and a 


virial radius of 2 pc, which result in a stellar density of about 
lOOstars/pc^ and a velocity dispersion of (w) ~ 2km/s. The 
lifetime of a star cluster in the Galactic disc can be estimated 
from (Lamers et al.|2005|: 


tci = 2.24 Myr 


[Ms] 


(4) 


For the star cluster in which the Sun was born this results 
in a lifetime of about 200 Myr. 

Our adopted cluster lifetime and the assumption of a 
constant density within this period are very approximate. 
However, our intention is to estimate the relative importance 
of experiencing an encounter in the parental star cluster or 
in the Milky Way Galaxy, for which this approach suffices. In 
order to estimate the importance of these assumptions we in¬ 
tegrated the cluster mass and radius evolution from the sim¬ 
ulations by [Portegies Zwart et al)] ( |2001[ ). They aimed their 
simulations at mimicking Pleiades, Praesepe and Hyades, 
which have comparable initial conditions as the Solar birth 
cluster. According to these simulations, clusters with such 
parameters survive for as long as a Gyr during which the 
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cluster mass drops linearly with time. In this period the 
cluster expands by abont a factor of 3, ronghly proportional 
to the sqnare root of time. With this slightly more elaborate 
estimate we argne that we overestimate the exposure of en¬ 
counters in the star cluster by about a factor of two. This 
would be consistent with adopting 100 Myr for the cluster 
exposure calculation in Figure[^ rather than 200 Myr. 

The distance to which these perturbations induced by 
close encounters penetrate into the Solar System, depends 
on the distance of closest approach gene, and the eccentricity 
of the encounter Bene, and can be expressed in (Kobayashi 
fc Ida| 2001 | |: 


u(l -|- e) ~ ( —) ^^genc ((1 + Af/M 0 )(l -|- Bene)) . (5) 

Recently Jilkova et al. (|2015| ) argue that the planetesi- 
mals Sedna ( Brown et al.||2004 l and 2012 VPii 3 (Trujillo & 


Sheppard 20141 were captured by a close encounter from 


a M ~ 1.8 Mq star that passed the Solar System with 
a closest approach of gene — 227 AU and relative veloc¬ 
ity V ~ 4.3km/s (which corresponds to an eccentricity 
Cenc = 2.6). According to Eq.[^ such an encounter would 
have perturbed the Solar System to a distance of about 
36 AU. In Figure[^ we present the distance to which such 
and encountering star perturbs the planetesimal around the 
Sun. For completeness we include the two objects Sedna and 
2012 VPii 3 to indicate how dramatically the encounter with 
star Q perturbed the Solar System. 



a [AU] 


Figure 2. The fragility of the Solar System in the Galactic field. 
The shades represent rienc: the probability density distribution 
for the number of encounters that affect the orbits in the Solar 
System to a given semi-major axis a and eccentricity e. Here 
we adopted a stellar density of n = 0.20 stars/pc®, a velocity 
dispersion of u = 30km/s and the typical mass of an encountering 
star of M = 0.5 Mq. The solid black curve and the two bullet 
points give the Neptune’s perturbing distance, and the orbital 
parameters for Sedna and 2012VPii3, see also Fig.^ The yellow 
curve gives the distance to which the Solar System was perturbed 
(by 1 % of its velocity at apohelion) due to the recent encounter 
with the 0.15 Mq binary star WISE J072003.20-084651.2, which 
grazed the Solar System, but clearly did not come in close enough 
to perturb the inner Oort cloud. 


2.3 The effect of external perturbations from the 
Galactic encounters 

Once the parental star cluster dissolves, the planetary sys¬ 
tem can only be perturbed by internal reorganization, and 
by the Galaxy. This latter can be subdivided in a global 
perturbation from the slowly varying Galactic tidal field, 
but the occasional close encounters with field stars are more 
important (see however |Fernandez|1997| ). 

The probability of an encounter with a Galactic field 
star is much smaller than of a close encounter in a star clus¬ 
ter, but the lower encounter rate is compensated in part by 
the longer time spent in the relatively low-density environ¬ 
ment of the Galaxy compared to the time spent in the star 
cluster. 

For estimating the effect of an encounter with a field 
star we cannot simply adopt Eq[^ because these are tuned 
for low-velocity (and low eccentricity) encounters, whereas 
Galactic encounters tend to occur with a much higher ve¬ 
locities. We therefore adopt the classic gravitatioal focussed 
cross section (Binney fc Tremaine||198T I 


1 + 


2G(m -b M) 


( 6 ) 


for calculating the encounter rate between the Solar System 
(m = IMq) and another Galactic star of mass M to a closest 
approach distance a. 

To provide an upper limit to the effect an encounter has 
on the orbital parameters of the planets or planetesimals we 
assume that the perturbed object and the closest approach 
are aligned. The heliocentric impulse gained by the object 
at distance r from the Sun is then given by (Rickman|1976 1 : 


Aw = 


2GM 


gene(gene w) 


(7) 


Here we assumed that the relative velocity, Wene = 30km/s, 
remains constant during the encounter and the mass of an 
encountering star M = 0.5 Mq. The perturbation is effective 
at apohelion (i.e., r = a(l-|-e)) and we estimate the distance 
to which the perturbation penetrates the Solar System at the 
point where the impulse gained by the object is comparable 
to its velocity at apohelion. This is a rather arbitrary choice, 
but suffices to indicate to which distance a passing field star 
may have affected planetesimals in the Solar System. 

In Fig.|^we present the number of encounters that per¬ 
turbed the Solar System by a passing Galactic disk star, 
after the parental star cluster has been dissolved. The prob¬ 
ability distribution is calculated using Eq.|^ and adopt¬ 
ing n = 0.20 stars/pc® and a velocity dispersion in the 
local standard of rest of w = 30km/s ( [Holmberg et ah' 
20071, and the time the Solar System spent in the Galaxy, 
foai = 4.3 Gyr. For comparison with Fig.[^ we include the 
curve for planetary perturbations (solid black curve) and the 
two objects Sedna and 2012 VPii 3 . 

Recently the Solar System had a close encounter with 
the M ~ 0.15 Mq binary star WISE J07200 3.20-084651.2 
(nicknamed “Scholz’s star” after its discoverer Schqlz|2014 ( 
at a distance of g = 0.251q o 7 pc (Mamajek et al.|2015 I. We 
estimate the perturbation by this encounter using the im¬ 
pulse approximation ( Heggie|1975 Rickman|1976 1 , in which 
the duration of the encounter is assumed to be much shorter 
that the period of the perturbed orbit (opposite to the adi¬ 
abatic regime in which Eq. is valid). The impulsive ap- 
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proximation is fulfilled for the encounter with Scholz’s star, 
which had a relative velocity of Venc = 83.2 km/s (Mamajek 


et al. 20151. We indicate this distance by the yellow curve 

Fig^ 

According to our analysis such an encounter (at gene ~ 
5.2-10"^ AU) should occur ~ 5000 times during the ~ 4.3 Gyr 
sojourn of the Solar System through the Galactic disk. Such 
encounter is therefore quite a likely event, which occurs 
roughly once every million years. The close approach of 
Scholz’s star occurred only 70,000 years ago ([Mamajek et al.| 


2015), which seems amazingly recent. If we naively divide 


the two time scales, such an encounter should already have 
happened ~ 60, 000 times. 

The e ffect o f this particular encounter has hardly per- 
(19271 cloud down to a distance of 10® AU 


Oort 


turbed the 

from the Sun. But an encounter with an equally low-mass 
star three orders of magnitude closer in would have affected 
the outer most planets. With the derived probability distri¬ 
bution Eq.|^ using Eq.j^ for the cross section and Eq.[^ to 
estimate the distance to which such an encounter affects the 
solar system. Such an encounter would be very unlikely to 
happen over the lifetime of the Solar System. 

Considering the analysis of jjilkova et al.j ( |2015| ) the en¬ 
counter that introduced the Sednitos (a family of planetes- 
imals with orbits similar to the 2003VBi2 Sedna) into the 
Solar System occurred in the parental cluster, and since then 
no other stellar encounter has perturbed Edgewordt-Kuiper 
( Edgeworth] 1943 Kuiper 1951) belt. This picture is con¬ 
sistent with the presence of a Parking zone (see § |3.1| ) in 
the Solar System between about 100 AU and 1000 AU; en¬ 
counters that affect the Solar System to a distance within 
1000 AU are extremely rare. 

The Parking zone in the Solar System is populated by 
the Sednitos ( Jilkova et al.|2015 I; planetesimals that share a 
common argument of pericenter, inclination and perihelion 
distance. We argue that the Parking zone in the Solar Sys¬ 
tem extends from about 100 AU to the about ~ 1000 AU, 
near the outer boundary of the inner Oort cloud. The distri¬ 
bution of orbital parameters of planetesimals discovered in 
this regime bear the information of the last strong encounter 
the Sun experienced from the time when it was part of its 
parental cluster. 


3 PERTURBING PLANETARY SYSTEMS IN 
GENERAL 

In Pig.[^ we present a generalized view of the fragility of 
planetary systems. The solid black curve gives the range 
to which possible massive planets perturb the inner parts 
of the planetary system. Violent planet scattering can cause 
massive planets to migrate further outwards, where they can 
perturb the local planetesimals ( Chatterjee et aL]|2008 l. In 
that case, the solid curve will shift to the right. 

The parental star cluster perturbs its planetary systems, 
but the further evolution of the planetary system determines 
to what range such a perturbation is preserved over time. 
For a globular cluster the stellar density is generally higher 
than for an open cluster, and the probability of spending a 
prolonged period in a globular cluster is also larger. As a 
result the range to which the planetary system is perturbed 
when born in a globular cluster is much closer to the star 


than for an open cluster. In fact, from the schematic pic¬ 
ture (Fig.|^ the influence of random encounters while the 
member of a globular clusters penetrates all the way to the 
inner planets. As a consequence, planetary systems in glob¬ 
ular clusters, or other massive dense star clusters, are likely 
to be perturbed by internal as well as external effects. This 
may explain, in part, the lack of observed planets in globular 
clusters ( Weldrake et al.|2007 |. 

If born in a low density cluster with a relatively short 
lifetime, the range to which random encounters penetrate 
into the planetary system hardly reaches the influence range 
of the giant planets. Once the planetary system escapes the 
star cluster, the perturbations from the Galaxy start to af¬ 
fect the orbits of the objects. This influence prolongs for the 
remainder of the main-sequence lifetime of the parent star, 
after which stellar evolution starts to play a major role in 
the redistribution of the orbits. 


3.1 The Frozen and Parking zones 

We define the Parking zone as a range in semi-major axis 
and eccentricity in which the orbits of objects have only 
been affected by encounters in the parental star cluster, 
and not by the local planets or by the Galaxy. 

Objects that orbit in the Parking zone have therefore not 
been affected by the planets and remain unaffected by close 
Galactic encounters. The Parking zone is likely to shrink 
with time, and young stars tend to have a more extended 
Parking zone than older stars, due to the less prolonged 
exposure to Galactic encounters. 

Objects with orbital parameters in the Parking zone 
preserve information about the last event that affected their 
orbits in the planetary system. Once in the Parking zone, 
orbital parameters are unlikely to be affected either by the 
planets, because they only affect orbits closer to the star, 
or by random Galactic encounters, because they tend to af¬ 
fect the outer most regions of the planetary system. Objects 
found in the Parking zone can therefore be used as tracers to 
reconstruct the event that introduced them in their current 
orbits. 

To the left of the Parking zone, and to the right of the 
range to which planets perturb the orbits we recognize the 
Frozen zone. 

We define the Frozen zone as a range in semi-major axis 
and eccentricity in which the orbits of objects have not 
been affected by the local planets and not by any encoun¬ 
ters, in the parental star cluster or the Galaxy. 

In this zone minor bodies remain unaffected by either in¬ 
ternal influences or external perturbations. Planetesimals 
found in this regime will preserve information about the for¬ 
mation of the planetary system. 

In the Solar System, the Frozen zone is probably very 
small or completely absent (see also Fig.[^. But other plan¬ 
etary system may have a populated Frozen zone, which can 
be used to study their origin. 
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Figure 3. The fragility of a planetary system under the influence 
of interaction with stars in its parental cluster and a galaxy. The 
horizontal axis has no scale, because the range to which the effect 
penetrates depends on the parameters of the planetary system 
and the environment. The black, yellow and blue curves give the 
ranges to which various perturbations affect the orbits of the mi¬ 
nor bodies. Specific example of the Solar System are presented in 
Figs.[^and|^ We identify the Frozen zone and the Parking zone. 
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